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ABSTRACT 

We present evidence that the high- velocity CivA1549 emission line gas of radio-loud quasars 
may originate in a disk-like configuration, in close proximity to the accretion disk often assumed 
to emit the low-ionization lines. For a sample of 36 radio-loud z «2 quasars we find the 20-30% 
peak width to show significant inverse correlations with the fractional radio core-flux density, R, 
the radio axis inclination indicator. Highly inclined systems have broader line wings, consistent 
with a high- velocity field perpendicular to the radio axis. By contrast, the narrow line-core shows 
no such relation with R, so the lowest velocity C iv-emitting gas has an inclination independent 
velocity field. We propose that this low-velocity gas is located at higher disk-altitudes than 
the high-velocity gas. A planar origin of the high-velocity C iv-emission is consistent with the 
current results and with an accretion disk-wind emitting the broad lines. A spherical distribution 
of randomly orbiting broad-line clouds and a polar high-ionization outflow are ruled out. 

Subject headings: galaxies: active — quasars: emission lines — quasars: geometry 



1. Introduction 

Little is currently known about the detailed 
geometry and structure of the broad-line region 
(BLR) in quasi-stellar objects (QSOs) due to the 
difficulty in spatially resolving the central regions 
even in the most nearby active galactic nuclei 
(AGN), the limited sensitivity of photo-ionization 
processes to the detailed BLR geometry, and to 
the complex nature of reverberation transfer func- 
tions (Peterson 1997; Netzer & Peterson 1997). 
For example, it is unclear whether the BLR gas has 
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a spherical or disk-like distribution, and whether 
the BLR consists of discrete clouds or is part 
of a continuous gas distribution, such as an ac- 
cretion disk wind [e. g., Boroson & Green 1992 
(BG92); Murray & Chiang 1997]. An improved 
knowledge of the geometry and structure of the 
centrally emitting and absorbing gas regions and 
their relation to one another and to the continuum 
source(s) may help us to better understand the 
observed continuum-line and line-line correlations 
(e. g., BG92; Brotherton 1996; Wilkes et al. 1999; 
Wills et al. 1999), and furthermore shed light on 
the interrelation of individual AGN subclasses. 

The best constraints on the geometry and struc- 
ture have so far been obtained from reverbera- 
tion mapping (see the excellent reviews by Pe- 
terson 1993; Netzer & Peterson 1997) of low- 
luminosity AGN showing (1) the BLR is rather 
compact (the outer radius, R ou ter is smaller than 
0.3 pc w 1 /larcsec), but is geometrically thick with 
Router ~ 10 Rinncr, (2) a certain degree of ioniza- 
tion stratification is present in the BLR gas, and 
(3) the BLR size and the (bolometric) luminosity 
of the active nucleus are correlated (e. g., Kaspi 
et al. 2000). Some clues to the relative position 
of the continuum source, and the scattering, emit- 



ting, and absorbing components, respectively, can 
be obtained from polarization studies of broad- 
absorption line QSOs (e. g., Ogle 1997) which ar- 
gue that they are highly inclined "radio-quiet" 
QSOs. Additional and independent clues to the 
BLR gas kinematics may be obtained by studying 
the widths of the broad emission lines themselves. 
The gravitational potential of the central engine 
is thought to dominate the kinematics of the BLR 
gas and the line width (FWHM) to estimate the 
gravitational potential at the location of the gas, 
(e. g., Peterson & Wandel 1999). If the motion of 
the BLR gas has a preferred direction, as opposed 
to disorganized virial motions around the central 
black hole, we would expect the projection of this 
velocity field on our sight line to vary as the source 
axis inclination varies in the AGN population. 

For a sample of broad-line radio sources Wills 
& Browne (1986) found FWHM(H (3) to relate to 
the radio core- fraction, log R, an estimator of the 
radio axis inclination to our line of sight, such 
that face-on sources lack very broad lines. They 
conclude that H/3 is emitted by a disk configu- 
ration perpendicular to the radio axis, which of- 
fers support to the two-component model for the 
high and low ionization lines (Collin-Souffrin et al. 
1988). The low-ionization, very high-density en- 
vironment of the accretion disk fits the require- 
ments of the physical conditions in the Balmer 
line emitting gas. The low level of long-term 
variability in the broad H/3 wings (e. g., Maoz 
et al. 1994; Peterson et al. 2000) seems to sup- 
port the disk interpretation [e. g., Corbin & Smith 
2000; but see Brothcrton (1996) for counter ar- 
guments]. Also, Zheng (1992) and Corbin (1997) 
model the high- velocity BLR emission finding con- 
sistency with Keplerian motion in low-z sources. 
Collin-Souffrin et al. suggested a polar outflow ori- 
gin for the high- ionization lines (HILs), thereby 
explaining their observed peak blueshift relative to 
the low-ionization lines (LILs; e. g., Gaskell 1982; 
Wilkes 1984; Tytler & Fan 1992). However, in the 
absence of a definitive model a spherical distribu- 
tion of dense, discrete HIL-emitting clouds is often 
adopted in BLR studies to explain trends such as 
the inverse correlation between [O in] A 5007 and 
the Fen emission (BG92). 

In order to explore these possible models we 
studied the inclination dependence of the widths 
of CivA1549, representative of the HILs, relative 



to those reported for H (3 in a sample of radio- 
powerful QSOs (RLQs). We find interesting con- 
straints applicable to the HIL geometry, which we 
present below. This project is part of an extended 
study of the ultraviolet (UV) broad emission lines 
in moderate- z QSOs with a range of radio power. 

2. Data Analysis 

The details of the sample and the optical data 
of the extended study are described by Vester- 
gaard et al. (in preparation, hereafter paper 1). 
We summarize the main points below. The ra- 
dio data are described by Barthel et al. (1988), 
Lonsdale, Barthel & Miley (1993) and Barthel, 
Vestergaard & Lonsdale (2000). The current sam- 
ple consists of 25 lobe-dominated (LDQs; R5 < 
0.5; R5 is defined below) and 11 core-dominated 
(CDQs; R 5 > 0.5) RLQs, chosen from the larger 
sample of z w 2 RLQs, originally selected from 
Barthel, Tytler & Thomson (1990), and the 3C 
and 4C catalogues. To ensure that we study the 
inclination dependence in the line widths, and 
not, say, a dependence on intrinsic source bright- 
ness, the CDQs and LDQs were selected such that 
their distributions of extended 5 GHz radio lumi- 
nosity, L ext , are similar. A Kolmogorov-Smirnov 
test confirms this to a confidence level of 97.7%. 
The L cxt -ranges and the mean and median values 
also match to within the standard error, while the 
My-ranges match to within the expected uncer- 
tainty of -0.5 m (Table 1). Three of the LDQs are 
compact (<15kpc) but their i?5-values could be 
measured from our high-resolution radio images. 
Highly variable RLQs are excluded. For the RLQs 
not presented by Barthel et al. (1990) we obtained 
MMT Blue Channel spectra during 1994 October 
- 1996 November. The RLQs were observed at air- 
masses below 1.4 and with a 1" slit placed at the 
parallactic angle. The Barthel et al. data were not 
well flux-calibrated and so were recalibrated using 
short MMT exposures. 

Spurious noise spikes and narrow absorption 
lines superposed on the emission lines were cor- 
rected by generating smooth fits to the profiles. 
The multiple Gaussian fitting (Laor et al. 1994; 
paper 1) that was applied to the Civ-Hen line 
complex also allows a deblending of the line emis- 
sion contributions, especially in the red wing of 
C iv, which in a few objects shows signs of weak 
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Fen emission. This is not expected to affect 
the current results adversely. The limitations im- 
posed by the Fe n emission will be discussed in pa- 
per 1. The profile parameters were measured on 
the smooth C iv fits to ensure that they are mea- 
sured consistently, objectively, and with minimum 
uncertainties. The full line widths at fractional 
profile heights, Y (= 10%, 20%, 90%) were de- 
termined (FWYM) , along with the interpercentile 
velocity widths, IPV(Y%) (Whittle 1985). Reli- 
able error estimates of each line parameter, based 
on continuum and flux uncertainties, were deter- 
mined in a similar fashion as outlined by Whittle 
(1985; see paper 1 for details). 

The radio core-fraction is defined here as R 5 
— Scorc/Stotai, where S„^f Q are the observed 
fluxes at 5 GHz measured from our VLA im- 
ages; the spectral index, a, is defined as a^oooMHz 
(or cksooomhI if thc former is unavailable). A K- 
correction 6 was applied to R 5 to permit a com- 
parison in the rest-frame: i? 5 = 5 , C orc/[(5'corc + 
S e xt.)(l + z) {aext -"—)], where a corc =0 and 
attended = -0.8 (e. g., Kellerman & Pauliny- 
Toth 1981). Our definition of R 5 is less sensitive 
to variations in thc density of the environment 
of the radio source than e. g., the R defined as 
Score/ S e xt., especially for face-on objects, where 
S'ext. tends to be relatively weak. We also con- 
sider the inclination measure, logi?v, defined as 
the ratio of the 5 GHz core flux density to the vi- 
sual flux density with applied K-correction (Wills 
& Brotherton 1995). It is also less sensitive to the 
density of the radio source environment. 

3. Trends with Source Inclination 

FWHM(Civ) shows no relation with logi? 5 
(the confidence level 7 [hcrafter CL] of the Spear- 
man rank correlation test is ~79%). However, 
FW20M(C iv), representative of the width of the 
broad C IV base, shows a strong inverse correlation 
with logi? 5 at the 99.7% CL (Fig. la). Support 
for the reality of this effect comes from the sim- 
ilar relations for FW10M, FW30M, and FW40M 
(CLs: 97.5%, 99.6%, and 98.6%). As an "area 
parameter", IPV(Y) is much less sensitive to the 
varying line-core strength which causes scatter in 

6 Ho = 50 km s _1 Mpc -1 , qo = is used throughout 
7 Confidence level = 1 — P, where P is the probability of no 
correlation 



the FWYM measures (paper 1). The fact that 
IPV(40-70%) show an inverse correlation with 
logi?5 at thc 97-98% CL, also supports the re- 
ality of thc FW20M trend. Thc IPV(Y<40%) is 
more sensitive to variations in the extreme wings 
than FWYM, evident from thc larger scatter and 
errors. This does not jeopardize the above results 
for IPV(40-70%) [paper 1]. Thc logi? v measure 
shows similar relations as log R5 though with rela- 
tively lower significance. This is most likely due to 
the fact that a zone of avoidance is more evident 
in these relations (e. g., Figure lb) though it is less 
apparent than that presented for H [3 by Wills & 
Browne (1986). Brotherton (1996), however, finds 
clear anti-correlations for FW25M and FWHM of 
H/3 (but none for FW75M). 

The Civ FW20-30M relations with core- 
dominance must primarily be driven by inclina- 
tion and not some other property because: (1) 
the L ex t distributions are similar for LDQs and 
CDQs. This is further confirmed by the at most 
weak FW20-30M - L cxt relations (CL <95%), (2) 
the lack of EW correlations with core-dominance 
(CLs <92%) indicates that the (luminosity nor- 
malized) line flux does not change with inclina- 
tion, (3) the lack of FW70-90M relations with 
core-dominance (CLs <37%) shows there is no 
weakening of the narrow line core in LDQs which 
could yield an apparent but unreal FW20-30M in- 
clination dependence, and (4) no correlation with 
inclination is found for My, L cx t, and redshift. 

The line width relations with inclination di- 
rectly show that the high-velocity field projected 
along the line of sight increases with source incli- 
nation. The velocity field of BLR gas (in discrete 
clouds or a continuous distribution) in organized 
orbits around the source axis is expected to behave 
this way. The similar H /3 and C IV relations im- 
ply that, to first order, the kincmatical structures 
of thc two high-velocity line-emitting regions, and 
by inference those of HILs and LILs, are similar. 
Item (2) above also implies that the line and con- 
tinuum fluxes are anisotropically emitted to about 
the same degree. Highly beamed sources are ex- 
cluded and thc UV continuum is also not expected 
to be significantly beamed (e. g., Wills & Broth- 
erton 1995). Item (2) confirms this. 

Though these results do not provide indepen- 
dent evidence for a disk-like geometry, they are 
consistent with the high-velocity HIL gas being 
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physically close to the accretion disk and with its 
predominant motion being in a plane parallel to 
the disk. However, the C iv emission cannot orig- 
inate in the low-ionization accretion disk-material 
itself, while the H/3 emission may (e. g., Collin- 
Souffrin et al. 1988). We therefore propose that 
the high- velocity HIL gas is planar and lies in close 
proximity to the disk, and with similar overall 
kinematics. In our ongoing work we will re-address 
this FW20-30M inclination dependence issue with 
larger and optimally constructed L oxt -samples of 
LDQs and CDQs. 

Our results argue against two possible BLR ge- 
ometries: (i) the velocity field of a spherical distri- 
bution of randomly orbiting clouds is not expected 
to vary with inclination, and (ii) a polar outflow of 
the HIL gas would show the opposite inclination 
dependence to that seen here for C iv. However, a 
planar HIL velocity-field and geometry is consis- 
tent with a disk-wind, where the broad lines are 
emitted at the base of dense, optically thick winds 
evaporating off the accretion disk (Murray & Chi- 
ang 1997, 1998). This model can successfully ex- 
plain a number of observational details. Bottorff 
ct al. (1997) and Konigl & Kartje (1994) present 
alternative disk- like BLR models. 

Item (3) above implies that the lowest- velocity 
BLR gas has a random velocity field relative to 
inclination. It may not be physically associated 
with the high-velocity gas located close to the ac- 
cretion disk. Low-density HIL gas of lower optical 
depth at higher disk-altitudes, similar to the pro- 
posed intermediate line region (ILR; e. g., Broth- 
erton et al. 1994), fits this requirement; the alti- 
tude is low enough to maintain the high ionization 
level. In the Murray & Chiang (1997) disk-wind 
model the narrow line-core is emitted by the wind 
at much larger disk-radii than the broad line-base. 
This predicts similar inclination dependences for 
the high and low velocity fields in contrast to our 
results and so requires a modification of the model, 
possibly along the lines suggested above. 

Brotherton (1996) finds very similar results for 
the inclination dependence of the high and low 
velocity H/3 emission in a sample of RLQs. The 
VBLR/ILR picture, which he advocates, is not ge- 
ometry specific, but is fully consistent with the 
results and interpretations presented here. 

A disk-like geometry of the HIL gas predicts: 
(1) other HILs (Lya, Nv, Siiv+Oiv], and Hen) 



should show similar dependences on source incli- 
nation; such a study is underway, (2) the low- 
velocity dispersion (< 2000 kms -1 ) of the nar- 
rowest line core of the HILs should be indepen- 
dent of inclination, assuming this gas is located 
at relatively high disk-altitudes, and (3) the Civ 
and H /3 profile parameters should behave similarly 
with varying source inclination, excluding the C iv 
narrow line core and the contribution to H (3 from 
the narrow line region. No detailed study exists of 
high-rcdshift QSOs; the FWHMs of Civ and H/3 
correlate somewhat at low redshift (e. g., Wilkes 
et al. 1999). We are currently studying this issue 
for low-rcdshift QSOs, based in part on Hubble 
Space Telescope archival data. 

4. Summary 

Despite years of study little is known for certain 
about the geometry of the BLR of AGN beyond its 
compact size (< light-weeks). The present results 
indicate that the high-velocity HIL gas is located 
in a disk-like configuration close to the accretion 
disk, long thought to emit the LILs, and with sim- 
ilar kinematics. Our results argue against a spher- 
ical cloud distribution with random cloud veloci- 
ties and a high-ionization polar outflow. If con- 
firmed, the implications are intriguing; they may 
help to constrain interpretations of BLR data and 
to suggest new studies to shed light on the detailed 
structure and kinematics of the AGN BLR. 
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Table 1 
Subsample Match Statistics. 
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L ext (5 GHz) 
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Mean 
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Range 


LDQ 
CDQ 


26 
11 


27.75 
27.78 


27.75 
27.82 


26.84 26.89 - 28.13 
27.12 27.10 - 28.11 


-26.5 - -28.9 
-27.0 - -29.2 
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Fig. 1— The FW20M relations with (a) logi? 5 
and (b) log-Ry- Core-dominated RLQs (CDQs) 
have solid symbols. Spearman rank correlation 
tests show confidence levels of (a) 99.7% and (b) 
96%. The lower confidence levels for logi?v is 
likely due to its relation being a zone of avoid- 
ance of very broad lines in CDQs and not a clear 
inverse correlation. 



